JournaloflnsectScience.Vol.il | Article 173 Chandra et al. 




Genetic variability in populations of the southern chinch bug, 
Blissus insularis, assessed using AFLP analysis 

Ambika Chandra' 3 *, James A. Reinert" 3 , Jonathan LaMantia lc , J. Blake Pond l,3d , David R. 
Huff 26 

'Texas AgriLife Research, Texas A&M System, 17360 Coit Road, Dallas TX 75252-6599 
2 Crop and Soil Sciences Department, Pennsylvania State University, University Park, PA 16802 
3 Current address, UT Southwestern Medical Center, 5323 Harry Hines Blvd. Dallas, TX 75390 

Abstract 

Southern chinch bug, Blissus insularis Barber (Heteroptera: Blissidae), is the most destructive 
insect pest of St. Augustine grass, Stenotaphrum secundatum Waltz (Kuntze), in the southern 
United States. The present study is focused on assessing genetic variability in five populations of 
B. insularis collected from Texas and Florida where St. Augustine grass is widely grown. The 
amplified fragment length polymorphism technique was used to DNA fingerprint individuals 
from each population (a total of 46 individuals) using five primer combinations {EcoRllMSeY). 
Analysis of molecular variance results show no evidence to support significant genetic variability 
among Texas and Florida populations of B. insularis. Nearly all genetic variation was found to 
reside within populations (95%), with only approximately 3% residing among populations 
between regions. Low Gst values obtained from POPGENE and low Fst values obtained from 
the analysis of molecular variance both support the conclusion for high levels of gene flow 
resulting from interbreeding and/or migratory events among the populations. A Mantel test of 
Euclidean squared distances showed no correlation between the genetic distance and geographic 
distance matrices of tested populations of B. insularis. The results of the present study suggests 
that gene flow is occurring among populations of B. insularis and, therefore, breeders need to be 
aware that this resistance will most likely not remain localized, and it has the potential to spread 
as a result of migratory events. 
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Introduction 



Chinch bugs (Hemiptera: Blissidae: Blissus), 
B. insularis, are common pests of agronomic 
crops and turfgrasses. The four most 
economically important species of chinch 
bugs prevalent in the United States are the 
southern chinch bug, B. insularis Barber; the 
common chinch bugs, B. leucopterus 
leucopterus; the hairy chinch bug, B. I. hirtus; 
and the western chinch bug, B. occiduus. 
Although germplasms resistant to individual 
species of chinch bugs have been identified 
(Reinert et al. 2004), a comprehensive study 
conducted by Anderson et al. (2006) 
documented the importance of developing 
grass varieties that exhibit resistance to 
multiple species of chinch bugs for wider 
adaptability, reducing the need for pesticides 
and promoting the adoption of integrated pest 
management strategies. 

All four economically important species of 
chinch bugs have been reported to have 
extensive and overlapping host range and 
geographic distributions (Anderson et al. 
2006; Pierson et al. 2007; Reinert unpublished 
data). These overlapping ranges warrant the 
need to implement molecular tools, such as 
DNA markers, for species identification and 
to assess genetic structures of the genus 
Blissus. The amplified fragment length 
polymorphism (AFLP) marker system has 
been successfully used in several population 
genetics studies of insects (Reineke et al. 
1999; Gaviria et al. 2006; Clark et al. 2007; 
Kazachkova et al. 2007; Krumm et al. 2008), 
including chinch bugs. Pierson et al. (2007) 
investigated the genetic variation between two 
species of chinch bugs, namely B. occiduus 
and B. I. leucopterus, using AFLP markers 
and showed that the majority (63%) of 
variability is present within species and 
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approximately 37% variability is present 
between species. To our knowledge, the 
genetic structure of populations within and 
among species of chinch bugs has yet to be 
reported. 

In the present study, the level of genetic 
variability in B. insularis, which is considered 
a major pest of St. Augustine grass, 
Stenotaphrum secundatum Waltz (Kuntze), 
throughout the southern United States was 
investigated. B. insularis has also been 
reported on other warm-season turfgrasses 
such as bermuda grass, Cynadon dactylon; 
bahia grass, Paspalum notatum; centipede 
grass, Eremochloa ophiuroides; and zoysia 
grass, Zoysia spp. (Kerr 1966). 'Floratam', a 
cultivar of St Augustine grass, was 
extensively planted in the southern United 
States following its release in 1973 (Horn et 
al. 1973) primarily because it exhibited good 
resistance to B. insularis. Resistance of 
'Floratam' to B. insularis is categorized as 
antibiosis because of high chinch bug 
mortality (80% to >90% mortality within a 7- 
day feeding period) and low oviposition rates 
(Reinert and Dudeck 1974). However, after 14 
years of commercial production and planting, 
the resistance of Floratam to B. insularis 
began to breakdown (Busey and Center 1987). 
Since 1987, turfgrass breeders have been 
striving to develop cultivar(s) with resistance 
to newer biotypes of B. insularis. The 
cultivars, 'FX-10' and 'Captiva', were 
developed and released in Florida in 1993 and 
2008, respectively (Busey 1993; Nagata and 
Cherry 2003), and provided resistance to 
Florida populations of B. insularis (>90% 
mortality within 14 days) that had overcome 
the resistance expressed in Floratam (Busey 
and Center 1987; Nagata and Cherry 2003). 
Preliminary studies show that these three 
cultivars of S. secundatum (Floratam, FX-10, 



Journal of Insect Science | www.insectscience.org 



2 



Journal of Insect Science.Vol. I I | Article 173 

and Captiva) are susceptible to populations of 
B. insularis collected from south Texas with 
<20% of the chinch bugs killed by these 
cultivars within a 7-day feeding period in 
laboratory studies (Reinert et al., manuscript 
in preparation). Based on these herbivory 
studies, there appears to be differences among 
populations of B. insularis. 

The objective of the present study was to 
assess the genetic variability between different 
populations of B. insularis using AFLP 
marker system. AFLP was chosen for this 
study because it does not require any prior 
knowledge of the genome (Vos et al. 1995) 
and generates a large number of polymorphic 
loci which are advantageous for studying 
genetic diversity and population structures 
(Salvato et al. 2002; Campbell et al. 2003; 
Blanc et al. 2006). 

Materials and Methods 

Sample collection and DNA extraction 

Populations of B. insularis were vacuum 
collected using a modified Echo leaf vacuum 
(www. echo-usa.com) from S. secundatum 
lawns growing in three geographically 
different locations in Texas and two locations 
in Florida (Table 1). In the laboratory, a 
mechanical aspirator was used to separate the 
chinch bugs from plant material and debris. 
Collected chinch bugs are maintained in 95% 
ETOH and stored at -20° C to avoid DNA 
degradation. Twenty individuals were chosen 
randomly from each population. For DNA 
isolation a modified cetyltrimethylammonium 
bromide method was used (Black and Duteau 
1997). Quantification of extracted DNA was 
performed to a known concentration of X 
DNA by loading 1.0 pi of extracted DNA 
mixed with 1.0 pi of tracking dye to a 1% 
agarose gel, which was electrophoresed for 30 
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min. Ethidium bromide stained gels were 
visualized using the Gel Doc System. 

Amplified Polymorphic Length 
Polymorphism (AFLP) 

The AFLP procedure was performed in three 
steps: 1) DNA template preparation, 2) DNA 
template pre-amplification, and 3) AFLP 
selective amplification. 



Table 1. Locations 
Augustine grass. 


of collected B. insularis populations from St. 


Population II) 


Location 


Latitude 


Lon^itudc 


Nt 
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32°22.2659 - 
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95°39.9N66' 


29°72.3357' 
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Highland Park 


96°8 1.3390' 


32 S3. 23 16' 
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s: 163760' 


29°30.I435" 
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FL 




Fl. Lauderdale 


S0°25.3443' 


26° 11. 5302' 


II 


|N = number of bugs analyzed per 
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DNA template preparation. A total of 7 ul 
of genomic DNA (1.0 pg) was incubated with 
EcoRl and Msel restriction endonucleases for 
2.5 h at 37° C followed by 15 min at 70° C in 
a total volume of 12.5 pi of restriction 
digestion mixture per reaction, which 
contained 0.0625 pi of EcoRl (100,000 U/ml), 
0.125 pi of Msel (50000 U/ml), 1.25 pi of 
10X One-phor-all, 0.125 pi of BSA (10 
mg/mL), and 3.938 pi of autoclaved nanopure 
water to make up the final volume. A total of 
5 pi of adapter ligation mixture; consisting of 
0.5 pi of 10X T4 DNA ligase buffer, 0.15 pi 
of T4 DNA ligase (2,000,000 U/ml), 0.5 pi of 
EcoRl adapter mix (5 pM), 0.5 pi of Msel 
adapter mix (50 pM) (Table 2), and 3.35 pi 
autoclaved nanopure water per reaction; was 
then added into each restriction digestion 
mixture and incubated for 8-10 h at 25° C. 
After incubation, the samples were diluted 
with 135 pi of TE buffer and stored in a -20° 
C freezer. All restriction enzymes, buffers, 
and ligation components were supplied by 
New England BioLabs supplies 
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(www.neb.com). Adaptors were synthesized 
by Integrated DNA Technologies 
(www.idtdna.com) 



Table 2. List of adaptors and selective pr 


mers used to 


characterize AFLP band patterns in five populations of 8. insularis. 




\he\ site 


rYruRI site 


Adapter* 


5' -GAOGATGAGTC'CTGA- 3' 


5' -I TCGTAGACTGIGTACC - 3' 




5-1 U H A< it iAl I < AT -.V 


f \ \1 IT .liiAl l it 1 At - V 


Prc-sclectivc 






primer 


5' -G ATGAQTCCTOAGTAA- V 


y -GACTGCGTACCAATTC-V 


Selective primer 


5-Msei-CrC- y 


5' *6FAM-EcoR 1 - A AG- 3' 




5' -Mscl-CTT-3 - 


?■ •<■! \\I 1 ei.Kl \ \ti '■' 




: ; -\1.e -( \ \- \ 


5 '(,1 oiK |. Ml . V 




f -M»el A AA- .V 


>■ •(.[ eciRi- a(k • y 




y-Msel-CAT-.v 


5' *6FAM-P.coR 1 - AGO- 3' 


* labeled primer 



Pre-amplification of DNA template. A 22 ul 

master mix containing 2.5 (0,1 10X Taq buffer, 
0.1 ul of Taq polymerase (5 U/ul), 1.5 ul of 
MgCl 2 (25 mM), 0.3 ul of dNTP mix (10 
mM), 0.5 ul of EcoKl-C (10 uM), 0.5 ul of 
Msel-A (10 uM), and 16.6 ul of autoclaved 
nanopure water was added to 3 ul of diluted 
ligation product obtained from the previous 
step. Each reaction tube was then placed in a 
thermal cycler and amplified using 30 PCR 
cycles of 72° C for 2 min, 94° C for 30 s, 56° 
C for 30 s, and 72° C for 2 min followed by 
60° C for 10 min. The pre-selective PCR 
product was diluted with 1 67 ul of TE buffer 
and stored at -20° C or immediately used in 
selective PCR. The oligonucleotide primers 
used in the pre-amplification of the DNA 
template are complementary to the 
adapter/restriction site with Msel primer 
containing one selective nucleotide and EcoRl 
primer containing no selective nucleotide 
(Table 2). All PCR components were supplied 
from Applied Biosystems/Life Technologies 
( www. appliedbiosy stems . com ) . Primers were 
synthesized by Integrated DNA Technologies. 

AFLP selective amplification. A total of 7 ul 
of PCR master mix was made consisting of 
1.25 ul of 10X Taq buffer, 0.1 ul of Taq 
polymerase (5 U/ul), 1 ul of MgCl 2 (25 mM), 



0.15 ul of dNTP mix (10 mM), 1 ul of BSA 
(10 mg/ml), and 3.5 ul of autoclaved nanopure 
water per reaction. A selective primer master 
mix was made for each of the five primer 
combination (Table 2) used consisting of 2.5 
ul of Msel (100 uM) selective primer and 0.5 
ul of 6FAM labeled EcoRl (100 uM) selective 
primer. The PCR master mix and 3 ul of the 
primer master mix were then mixed with 2.5 
ul of pre-selective PCR products in PCR 
tubes. The tubes were placed in a thermal 
cycler and amplified using following 
parameters: 94° C for 2 min, 1 cycle of 94° C 
for 30 s, 65° C for 30 s (lowering the 
annealing temperature by 0.7° C per cycle), 
72° C for 2 min; then 23 cycles of 94° C for 
30 s, 56° C for 30 s, and 72° C for 2 min, and 
finally 72° C for 10 min before holding at 4° 
C. After selective PCR amplification was 
complete, the selective PCR products were 
stored at - 20° C or were immediately 
prepared for AFLP analysis. In preparation of 
AFLP analysis, 0.5 ul of GeneScan - 500 
ROX Size Standard (ABI), 8.5 ul of Hi-Di 
Formamide (Applied Biosystems), and 1 ul of 
selective PCR product were mixed and placed 
in 96-well plates. The plates were analyzed 
using the ABI3130 capillary system. 

Data analysis 

A total of 501 AFLP markers were scored for 
each individual bug. For each individual bug, 
AFLP markers were scored as Is when 
present and 0s when absent. Data was scored 
using Genographer version 2.1 using file type 
- Gene Mapper Filter and Sizing Algorithm - 
Cubic Spline Interpolator. Bands with a peak 
height 300 or greater and 70 bp or larger were 
scored. All band sizes were rounded down to 
the nearest whole basepair and placed in a 
basepair length bin. Bins with multiple bands 
were considered stutters and only counted 
once. All bins with bands were converted to 1 
and all empty bins were converted to 0. 
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Gene diversity and gene flow. The resulting 
data was analyzed using POPGENE version 
1.32 (Yeh and Boyle 1997) to assess the 
percent polymorphism, gene diversity 
estimates of total populations (Ht) based on 
Nei (1973), genetic differentiation between 
populations (Gst), and gene flow estimates 
(Nm). Assuming Hardy- Weinberg equili- 
brium, the analyses were conducted at three 
levels: individuals, regions, and whole 
population using dominant marker data. H T 
values were estimated for each individual 
population and Gst values were estimated for 
populations within region and all five 
populations. Gst is expressed as H T - H s / H T 
where H T is gene diversity of the total 
population and H s is gene diversity of the 
single population. Gene flow was estimated 
from Gst expressed as Nm = 0.5 (1- Gst / Gst) 
(McDermott and McDonald 1993). 

Analysis of molecular variance. Analysis of 
molecular variance (AMOVA) (Excoffier et 
al. 1992; Huff et al. 1993) was conducted 
using Arlequin 2.0 (Schneider et al. 2000) to 
asses the genetic structure and genetic 
variability present within and among 
populations. Total variation of the AFLP data 
set was partitioned into three components, 
namely: among regions, among populations 
within a region, and within populations. 
AMOVA uses a permutational testing 
procedure to determine the significance of 
these three variance components. Pairwise 
comparisons were conducted to test genetic 
divergence among populations (Fst). 

Phylogenetic analysis. A phylogenetic 
analysis was conducted using PAUP 4.0 Beta 
version (Swofford 1998). A neighbor joining 
tree was generated using PAUP based on the 
Euclidian squared pairwise genetic distance of 
the 46 individual bugs belonging to five 
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different regions obtained using Arlequin. 
MEGA version 3.1 (Kumar et al. 2004) was 
used to view and label trees generated from 
PAUP 4.0. Bootstrap analysis was conducted 
to evaluate branch robustness of the resulting 
trees. Bootstrap support was estimated using 
1000 bootstrap replicates, and each replicate 
consisted of 10 heuristic searches, random 
addition sequences with branch swapping, and 
the multree option off. 

Principle components analysis (PCoA) 

Principle components analysis was conducted 
using NTSYSpc version 2.1 (Rohlf 2002) to 
construct a three-dimensional plot of the data 
for better visual representation. 

Mantel test. A Mantel test (Mantel 1967) 
with 1000 permutations was conducted using 
the MxCOMP feature of the NTSYSpc 
version 2.1 to test the correlation between 
geographic distance and the pairwise genetic 
distance between the five sampled 
populations. 

Results and Discussion 

Individual bugs from each population that 
produced banding patterns using all five 
primer pairs were used in the AFLP analysis. 
This resulted in seven individuals from the 
Longview, TX population; ten from the Rice 
University, Houston, TX population; nine 
from the Highland Park, TX population; nine 
from the Gainesville, FL population; and 11 
from the Ft. Lauderdale, FL population - 
making a total of 46 individuals analyzed in 
the present study (Table 1). The five primer 
combinations used to analyze the 46 
individuals from the five populations of B. 
insularis generated a total of 501 polymorphic 
loci with their fragment sizes ranging from 
70-495 bp. High levels of polymorphisms 
were observed within each population ranging 
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from 68.66% for the Longview, TX 
population to 85.43% for the Highland Park, 
TX population (Table 3), with an average of 
78.04%. 



Table 3. Descriptive statistical analysis of AFLP markers, across 
five populations of 8. insularis belonging to two regions, using 
POPGENE version 1.32. 


Population ID 


Polymorphic 
loci (%) 


Hr(SD) 




Nm 


Texas 








0.09 1 8 


4.9482 




LV 


6X66 


0.20 ± 0.18 








Rt 


7S.24 


0.22- 0.17 








HP 


85.43 


0.26 ±0.17 






Florida 








0.0493 


9.6352 




GV 


73.25 


0.19 ±0.16 








FL 


84.63 


0.22 ±0.1 6 


















All five populations 


100 


0.24 ± 0.02 


0.0261 


18.6478 


Ht = Nei's (1973) total gene diversity estimates, SD = Standard 
deviation, Gst = magnitude of gene differentiation among 
populations, Nm = Gene flow estimate. 



Genetic diversity, measured as population 
heterozygosity (H T ) values calculated from all 
sub-populations, was high and ranged from 
0.1898 for the Gainesville, FL population to 
0.2574 for the Highland Park, TX population 
with an average gene diversity of 0.21844 
(Table 3). Analysis across all five of the 
populations also revealed high genetic 
diversity value of 0.2428 among populations. 
Genetic variation within and among 
populations was measured as Gst- Gst values 
of 1.0 indicate that a majority of genetic 
variability resides between populations while 
lower Gst values (<0.5) indicate that a 
majority of the genetic variability resides 
within a given population. In the present 
study, Gst value among all sampled 
populations was low (0.0261) indicating a 
high degree of genetic variability within 
population (97.4%) and low variation among 
populations (2.61%) (Table 3). These results 
indicate that there is apparently little genetic 
differentiation among the sampled 
populations. The lack of genetic 
differentiation among populations is generally 
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considered the result of sufficient gene flow, 
typically in the form of migration, occurring 
across populations to counteract any effects of 
selection and/or genetic drift. 

Gene flow is commonly measured as Nm 
where N is the number of individuals in a 
population and M is the proportion of 
individuals in the population as a result of 
immigration. Nm values >1 indicate a high 
level of gene flow and that the effects of gene 
flow on population differentiation is greater 
than the effect of random genetic drift 
(Bossart and Prowell 1998). In the present 
study, Nm values were high in two tested 
regions, Texas and Florida, (Table 3) 
suggesting high gene flow from events such as 
interbreeding and migration. High gene flow 
values also suggest that populations among 
species will become genetically homogeneous 
in the absence of counteracting forces such as 
strong differential selection. 



Table 4. Analysis of molecular variance 


(AMOVA) for five 


populations of 8. insularis using Arlequin 


version 2.0. 








Sum of 


Variance 


Percentage 






Source of variation 


d.r. 


squares 


luiiiponclils 
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P value 


F Statistics 


Among groups 




I3&S34 


l.5<>3<>7 


o.s 


0.09091 


1 , , O.OITOi 


Among populations 














within groups 


3 


2°7.462 


2.50386 


3.10 


0.0X1293 


F M =0.03 166 


Within populations 


4 


■ . i'' •--:>., 


-r. s-_im 


°4.°2 


l) l« i 


1 o.nslMI 


Total 


4s 


'!".< 543 


su.tr:ii 








. . 



AMOVA showed that the difference among 
groups was not significant (P = 0.09091), but 
that the amounts of genetic variation among 
populations within groups and within 
populations were significant at P = 0.00293 
and 0.0000, respectively (Table 4). AMOVA 
showed that approximately 95% of the 
variation in the data set was from genotypic 
variation within populations. Only 3% of the 
variation was attributed to differences among 
populations within regions while the 
remaining 2% was due to the variation among 
regions (Table 4). A Mantel test (Figure 1) 
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showed that the genetic and geographic 
distances matrices were not significantly 
correlated (N = 10; r = 0.22021; P = 0.7682) 
indicating that there was no genetic isolation 
in regions of southern chinch bug populations 
tested in this study. Genetic divergence (F S t) 
values less than 0.2 indicates high gene flow 
and in the present study F S t for five sub- 
populations was 0.05080 indicating low 
differentiation among populations (Table 4). 
These results are consistent with those 
obtained from POPGENE (Table 3). This low 
level of differentiation between populations 
supports the possibility of gene flow among 
all the populations examined. 
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o 

H 




Genetic Distance 



Figure I. Mantel test (N = 10; r = 0.22021; P = 0.7682) of 
genetic distance versus geographic distance indicating lack of 
correlation. High quality figures are available online 



As expected, a phylogenetic neighbor-joining 
tree based on Euclidian-squared, pairwise 
genetic distance of 46 samples of B. insularis 
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Longview. TX(4) 
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Rice Universty. TX ( I ) 
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Figure 2. A phylogenetic neighbor joining tree based on 
euclidian squared pairwise genetic distance generated using 
PAUP 4.0 Beta. 46 8. insularis belonging to five different regions 
included in the tree. Underlined samples are from the 
Longview, TX population. Boostrap values are indicated above 
branch. Number in Parenthesis indicates the number for the 
bug sample analyzed in each population. High quality figures are 
available online 



(Figure 2) shows no clearly defined genetic 
structure among regions. However, the 
Longview, TX population (underlined in 
Figure 2) appeared to be a particularly unique 
population where individuals from within this 
population tend to cluster with individuals 
from the Florida populations. This scenario is 
more apparent when principle component 
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analysis was conducted and a three- 
dimensional Eigen plot (Figure 3) was 
graphed, which depicts mean pairwise genetic 
distance for the five populations of B. 
insularis. This situation leads us to speculate 
that the Longview, TX population is likely 




Figure 3. Principle coordinate analysis based on mean pairwise 
genetic distance for five different population of 8. insularis using 
NTSYS pc22. High quality figures are available online 



leading to the non-significant genetic 
variability among regions in our AMOVA 
analysis. Therefore, the next step was to 
conduct an AMOVA analysis without the 
inclusion of the Longview, TX population. 
However, when an AMOVA was conducted 
without the inclusion of the Longview, TX 
population, the difference among groups was 
again non-significant (P = 0.36950) and the 
difference in genetic variation among 
populations within groups and within 
populations were both significant at P = 
0.03226 and 0.0000, respectively. Without the 
inclusion of the Longview, TX population, the 
AMOVA showed that 93.4% of the variation 
in the data set was from genotypic variation 
within populations, 2.4% of the variation was 
attributed to differences among populations 
within regions, while the remaining 4.2% was 
due to the variation among regions. Therefore, 
an AMOVA either with or without the 
inclusion of the Longview, TX population 
gave similar results, in that, there was no 
significant genetic variability present among 
regions, and that the majority of the genetic 
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AFLP variability measured was found to be 
present within populations. 

Based on the results of the present study the 
resistance mechanisms of S. secundatum to B. 
insularis do not appear to be related to 
specific population differences in B. insularis, 
but rather it is more likely that such resistance 
results from the expression and regulation of 
specific genes from individual bugs. The 
present study also suggests that gene flow is 
occurring among populations of B. insularis 
and, therefore, breeders need to be aware that 
this resistance will most likely not remain 
localized and it has the potential to spread as a 
result of migratory events. In the future, a 
more detailed investigation might be 
performed by increasing the number of 
populations of B. insularis from the 
geographic distribution of S. secundatum; 
especially - a few more populations collected 
and analyzed from the Longview, TX region 
might be informative and useful in better 
understanding the population genetic structure 
of B. insularis. 

Acknowledgements 

This work was performed using the new 
faculty financial support provided by the 
Texas AgriLife Research. 

References 



Anderson WGT, Heng-Moss TM, Baxendale 
FP, Baird LM, Sarath G, Higley L. 2006. 
Chinch bug (Hemiptera: Blissidae) mouthpart 
morphology, probing frequencies, and 
locations on resistant and susceptible 
germplasm. Journal of Economic Entomology 
99:212-221. 

Black WC, Duteau S. 1997. RAPD-PCR and 
SSCP analysis for insect population genetic 
studies. In: Crampton J, Beard CB, Louis C, 



rnal of Insect Science | www.insectscience.org 



Journal of Insect Science.Vol. I I | Article 173 

Editors. The molecular biology of insect 
disease vectors: a methods manual, pp. 361- 
373. Chapman and Hall. 

Blanc MP, Lugon-Noulin N, Panighini C, 
Pijnenburg H, Rossi L. 2006. Structure of 
worldwide populations of Lasioderma 
seerricorne as revealed by amplified fragment 
length polymorphism profiles. Bulletin of 
Entomological Research 96: 111-116. 

Bossart JL, Prowell DP. 1998. Genetic 
estimates of population structure and gene 
flow: Limitations, lessons, and new directions. 
Trends in Ecology and Evolution 13: 171-212. 

Busey P. 1993. Registration of 'FX- 10' St. 
Augustine-grass. Crop Science 33: 214-215. 

Busey P, Center BJ. 1987. Southern chinch 
bug (Hemiptera: Heteroptera: Lygaeidae) 
overcomes resistance in St. Augustine grass. 
Journal of Economic Entomology 80: 608- 
611. 

Campbell D, Duchesne P, Bernatches L. 2003. 
AFLP utility for population assignment 
studies: analytical investigation and empirical 
comparison with microsatellite. Molecular 
Ecology 12: 1979-1991. 

Clark PL, Molina-Ochoa J, Martinelli S, 
Skoda SR, Isenhour DJ, Lee DJ, Krumm JT, 
Foster JE. 2007. Population variation of the 
fall armyworm, Spodoptera frugiperda, in the 
Western Hemisphere. Journal of Insect 
Science 7:05, Available online, 
http://insectscience.org/7-05 

Excoffier I, Smouse PE, Quattro JM. 1992. 
Analysis of molecular variance inferred from 
metric distances among DNA haplotypes: 
application to human mitochondrial DNA 
restriction data. Genetics 131: 479-491. 



Chandra et al. 

Gaviria DA, Aguilar E, Serrano HJ, Alegria 
AH. 2006. DNA fingerprinting using AFLP 
markers to search for markers associated with 
yield attributes in the silkworm, Bombyx mori. 
Journal of Insect Science 6:15, Available 
online, http://insectscience.Org/6.15 

Horn GC, Dudeck AE, Toler RW. 1973. 
'Floratam' St Augustine grass: A fast growing 
new variety of ornamental turf resistant to St. 
Augustine decline and chinch bugs. Florida 
Agricultural Experimentation Station Circular 
S-224. 

Huff DR, Peakall R, Smouse PE. 1993. RAPD 
variation within and among natural 
populations of outcrossing buffalograss. 
Theoretical and Applied Genetics 86: 927 
934. 

Kazachkova N, Meijer J, Ekbom B. 2007. 
Genetic diversity in pollen beetles (Meligethes 
aeneus) in Sweden: role in spatial, temporal 
and insecticidal resistance factor. Agricultural 
Forest Entomology 9: 259-269. 

Kerr SH. 1966. Biology of the lawn chinch 
bug, Blissus insularis. Florida Entomologist 
49: 9-18. 

Krumm JT, Hunt TE, Skoda SR, Hein GL, 
Lee DJ, Clark PL, Foster JE. 2008. Genetic 
variability of the European corn borer, 
Ostrinia nubilalis, suggests gene flow 
between populations in the Midwestern 
United States. Journal of Insect Science 8:72, 
Available online, http://insectscience.Org/8.72 

Kumar S, Tamura K, Nei M. 2004. 
MEGA3: Integrated software for 
Molecular Evolutionary Genetics 
Analysis and sequence alignment. 
Briefings in Bioinformatics 5: 150-163. 



Journal of Insect Science | www.insectscience.org 



9 



Journal of Insect Science.Vol. I I | Article 173 



Chandra et al. 



Mantel NA. 1967. The detection of disease 
clustering and a generalized regression 
approach. Cancer Research 27: 209-220. 

McDermott JM, McDonald BA. 1993. Gene 
flow in plant pathosystems. Annual Review of 
Phytopathology 31: 353-373. 

Nagata R, Cherry R. 2003. New source of 
southern chinch bug (Hemiptera: Lygaeidae) 
resistance in a diploid selection of St. 
Augustine grass. Journal of Entomological 
Science 38: 654-659. 

Nei M. 1973. Analysis of gene diversity in 
subdivided populations. Proceedings of the 
National Academy of Science USA 70: 3321- 
3323. 

Pierson LM, Serikawa R, Heng-Moss TM, 
Foster JE. 2007. An investigation of the 
genetic variation between Blissus occiduus 
Barber and Blissus leucopterus leucopterus 
(Say). Review of Undergraduate Research in 
Agricultural and Life Sciences 2(1): Article 2. 

Reineke A, Karlovsky P, Zebitz CPW. 1999. 
Amplified fragment length polymorphism 
analysis of different geographic populations of 
gypsy moth, Lymantrai dispar (Lepidoptera: 
Lymantriidae). Bulletin of Entomological 
Research 89: 79-88. 

Reinert JA, Dudeck AE. 1974. Southern 
chinch bug resistance in St-Augustine-grass. 
Journal of Economic Entomology 67: 275- 
277. 

Reinert JA, Engelke MC, Read JC. 2004. Host 
resistance to insects and mites, A review - A 
major IPM strategy in turfgrass culture. Acta 
Horticultura 661: 463-486. 



Rohlf FJ. 2002. NTSYSPC: Numerical 
Taxonomy Systems, ver. 2.1 IT. Exeter 
Publishing Ltd. 

Salvato P, Battisti A, Concato S, Masutti L, 
Patarnello T, Zane L. 2002. Genetic 
differentiation in winter pine processionary 
moth (Thaumetopoea pityocampa - 
wilkinsoni complex), inferred by AFLP and 
mitochondrial DNA markers. Molecular 
Ecology 11: 2435-2444. 

Schneider S, Roessli D, Excoffier L. 2000. 
Arlequin: software for population genetics 
data analysis. Version 2.000. Genetics and 
Biometry Laboratory, Department of 
Anthropology, University of Geneva. 

SwoffordDL. 1998. PAUP*. Phylogenetic 
Analysis Using Parsimony (*and Other 
Methods). Version 4. Sinauer Associates. 

Vos P, Hogers R, Bleeker M, Reijans M, Van 
De Lee T, Homes M, Frijters A, Pot J, 
Peleman J, Kuiper M, Zabeau M. 1995. 
AFLP: a new technique for DNA 
fingerprinting. Nucleic Acids Research 23: 
4407-4414. 

Yeh FC, Boyle TJB. 1997. Population genetic 
analysis of co-dominant and dominant marker 
and quantitative traits. Belgian Journal of 
Botany 129: 157. 



Journal of Insect Science | www.insectscience.org 



10 



